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In a previous paper (Eur. Phys. J. B 30, 239-251 (2002)) we have presented the main features 
and properties of a simple model which -in spite of its simplicity- describes quite accurately the 
qualitative behaviour of a quantum wire. The model was composed of N distinct deltas each one 
carrying a different coupling. We were able to diagonalize the Ifamiltonian in the periodic case and 
yield a complete and analytic description of the subsequent band structure. Furthermore the random 
case was also analyzed and we were able to describe Anderson localization and fractal structure of 
the conductance. In the present paper we go one step further and show how to introduce correlations 
among the sites of the wire. The presence of a correlated disorder manifests itself by altering the 
distribution of states and the localization of the electrons within the system 

PACS numbers: 03.65.-w: Quantum Mechanics, 71. 23. An: Theories and Models; LocaUzed States, 73.21.Hb: 
Quantum Wires 



Introduction 

In a previous papeni, hereafter referred to as I, the 
authors have developed a simple model describing the 
main features shown by a Quantum Wire, namely: Band 
Structure when the structure is periodically arranged, 
Anderson Localization in the random case and easy cal- 
culation of Lyapunov exponents, density of electronic 
states and electrical conductance endowed with an en- 
couraging fractal behaviour. We address the interested 
reader to reference I for the details of the model and the 
intrincacies of the computer calculations. 

The main drawback of the scheme thereby presented 
was the lack of a coherent description of the possible cor- 
relations which must be present in any realistic model 
of one dimensional conductance. It is the purpose of 
the present paper to fill this gap by presenting a sur- 
prisingly simple way to incorporate correlations entirely 
based upon the random nature of the model. This means 
that in addition to describe correlations with different 
parameters of the wire (spacing, number and amount of 
species, ... ) the model is flexible enough to be suitable 
for applying it to different correlation schemes such as 
the Tight-Binding model, the Anderson model and likely 
many others correlation prescriptions with sensible phys- 
ical content. 

The paper is organized as follows. In Section 1 we ana- 
lyze the meaning of the correlations in terms of probabili- 
ties of a given atom to be followed by other of a different 
character. After this analysis we are able to set some 
limits for the values of these probabilities. The next step 
is to incorporate those probabilities to the Functional 
Equation developed in I. How this functional equation is 
modified by the introduction of correlations in the way 
just described and how it can be used to calculate the 
correlation-corrected Lyapunov exponent is the subject 
of Section 2. In Section 3 some results are shown along 
with a discussion about the effects of the correlations to 
close finally with a section of Conclusions. 



I. CORRELATIONS 



Let us consider disordered quantum wires with short- 
range correlations which will be specified by the proba- 
bility for a certain atomic species to be followed by an- 
other type of atom in the chain sequence. In this model, 
the system will be characterized by the species concen- 
trations {ci}i=i_....„i where m is the number of different 
species, and the set {py }ij=i,...,m where pij means the 
probability for i to be followed by j (the prob. of finding 
j right after i). The correlations introduced in this way 
must be the consequence of the existence of an atomic in- 
teraction (or the effect of some physical parameters such 
as the size of the atoms) which might choose certain spa- 
tial sequences of the atoms modifying subtly the other- 
wise chain's purely random character. Thus this proce- 
dure seems a natural manner to account for correlations 
which can be present in nature or even those that can 
be produced inside a manufactured disordered quantum 
wire in a nonintentional way aside from the fact that of 
course one can always try to construct a wire exhibiting 
the desired correlations to recover the theoretical results. 

Let us clarify the meaning of pij. It is the ratio of 
the number of -ij- clusters and the number of i atoms 
(Ni). The probability pij is not the probability for j to 
be preceded by i, because the latter would be obtained 
dividing by the number of j atoms (Nj) instead of Ni. 
In the limit N ^ oo, the system will be determined by 
{ci} and {pij} and therefore the density of states will be 
given by a limiting distribution depending only on those 
quantities. And that distribution is of course the same 
whether we move from left to right or from right to left 
along the chain. So, the correlations we measure have to 
coincide in both directions. Then it is clear that for the 
infinite chain pij means the probability for an 'i ' species 
atom to be followed or preceded by a 'j ' species atom. 

Once we have understood the meaning of the correla- 
tions, let us derive the relations among them. What is 
the probability of finding at any position of the chain the 
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cluster -iji-? We can write this quantity as, 

Cj • Pij ■ Pji — >prob. of finding i ■ finding j after i ■ 

■ finding i after j 

Cj ■ — >prob. of finding j ■ finding i after j ■ 

■ finding i before j. 

Thus the equations for the correlations are, 
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< P^J < 1 



(la) 
(lb) 

(Ic) 



According to the equations, the correlations matrix (pij) 
is completely known from the above diagonal elements, 
"^^^"^^ where m is the number of different species. How- 
ever these elements are not completely independent due 
to ((Tct . because when one of the correlations is chosen, 
the maximum allowed value for some of the rest may be 
affected. This fact can be clearly seen in two simple ex- 
amples: 

2 species: the matrix is completely determined by pi2, 

ci 



Pii = 1 - Pi2 



P21 



C2 



-Pl2 



but it must be p2i < 1 P12 < — 

Cl 



P22 

C2 



1 -P2I 



Therefore pi2 < min < 1, 
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Cl 

3 species: the matrix is determined by Pi2,Pi3,P23- 
And we choose their values in that order. 

P12 < min < 1, — > which implies 

I Cl J 

<- ■ /1 

P13 < mm < 1 -pi2, — 

I Cl 

On the other hand it must be, 

P23 < 1 - P21 
C3 

P32 < 1 -P31 P23 < —(1 -P3l), 



C2 

therefore P23 < min < 1 -pi2, — -pi3 

L C2 C2 C2 

An expression for the general form of the maximum 
values can be obtained for arbitrary m. Choosing the 
above diagonal elements of the correlations matrix by 
rows. 
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FIG. 1: Correlation space for 2 species as a function of the 
concentration. The blue dashed line corresponds to the com- 
pletely random configurations. 



then 

Pij £ min 



( ^ i-i 3-1 ^ ^ i-i 

I * k=l k=i+l * k=l 




where the limits are given by the previous chosen corre- 
lations. 

Once the concentrations and P are fixed, the charac- 
terization of the quantum wire can be given in a compact 
form using the matrix 



Q = = (CiPy), 



(3) 



which is symmetric (Q = Q*) and it determines uniquely 
the wire and its correlations: 



2=1 



qij 

P^j = — 



Qij means the probability of finding the cluster -ij- (-ji-) 
at any position inside the wire. 

This model naturally includes the situation in which 
the disorder in the wire is completely random^, that is 
just defined by the values pij — Cj i ~ I, . . . ,m. Figure 
Ushows an example of the correlation space for 2 species. 



II. THE FUNCTIONAL EQUATION 

We will now obtain the correlation form of the func- 
tional equation, which enable us to calculate the density 
of states. In this process we are using some of the results 
and expressions obtained in reference I. 

The relations for the phase transmission inside the wire 
read. 



Tf'^ie) = arctan 2/ij(e) 



1 



Tr\e + n7T)^T-\e) 



tan 6'^ 
9e [0,7r) 



(4a) 
(4b) 
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where hj{e) = cos(e) + sin(£). To carry out the 

phase average we rearrange the terms depending on the 
species of the last atom we have passed, 



m. 



where Ni is the number of atoms of species i and the 
sumatory over ji means we arc summing over all the po- 
sitions j such that the atom at (j — 1) is of i-kind: 

{AO) = ci{M)i + ... + c,„{A9)ra = ^c^(A0)^. (6) 

7 

From now on the greek indices will refer to species while 
the latin indices mean the sites of the chain. In this case 
there exist correlations in the sequence of atoms, that 
is the probability for one species to be a neighbour of 
another type can be different for all of them. This implies 
that each species will give rise to a particular distribution 
function for the phase. 



-T 



d0 



{T,i0)-e}de (7) 



where Wj{9) is the distribution function of the phase at 
the position j generated by the species 7. The equations 
for the distributions ((B.14) in I) become here, 

WJ{9) = {T-\0)) ~ W^*_\ (I) + 1 (8a) 

W]{e + rn) = Wj{e) +r 61 G [0, tt) r e Z (8b) 

for all species and positions, where the superscript "(*)" 
means that we do not know the species of the atom pre- 
ceding the 7 atom, placed at position (j — 1); wc only 
know the probability for a certain species to be there. 
Let us introduce now the average over all possible se- 
quences with the given concentrations. 



and. 



(9) 



(10) 



The next step consists in taking the limit ^ cx) to 
approach the limiting distribution. 

We define Wy{0) = hm — V W^'iO) which is the 
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solution of: 

(11a) 

W^{e + rn) ^W^{e) + r 6'e[0,7r) r G Z. (lib) 



And finally we have to calculate, 



(AO) = ""-rP-rP I 

7,/3 



^'^^^mo)~e}de. (12) 



To carry out the integral we shall make use of the ex- 
istence of a value 9^ such that Tj{9q) = 9i V7, and 
proceeding by parts. 



® 



{A9) = {9, -9o)-Y, C7P7/3 / "^^ W^(0) 



d9 



dO 



Vc^ / W^{0)dO. (13) 



To evaluate the middle term we use the average of equa- 
tions lllal) which is 



1,13 



(14) 

With a simple change of variable, the use of (|14|l and the 
equations for the correlations we can write. 
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l+TT 



[T-\9)]de. 



(15) 



and going back to ((T^ . 

{A9) = {9,-9o)+7T-7TY,cpWp 



j W^{9)d9-J W^{9)d9y (16) 

It is not hard to see with the help of Ijllbp that the last 
term is 9q — 9i. Therefore, 



{A9)ie) 



Thus the density of states reads ((B.8) in I), 

dW, (f ) 



de 



(17) 



(18) 



which can be calculated numerically solving cqs. Ullaj) 
for several values of the energy. 



A. The Lyapunov exponent 

The quantity that provides information about the de- 
cree of spatial localization of the electronic states is the 
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Lyapunov exponent, which gives the exponential growth 
rate of the wave function, and thus its inverse is a mea- 
sure of the localization length of a state exponentially 
localized inside the wire. The Lyapunov coefficient can 
be written for the infinite wire as 



lim — log 

j=0 



N 



J + 1 



(19) 



where is the amplitude of the wave function at the 
jth atomic site. Using a mapping technique originally 
proposed in Ref[3 we are able to write A as a function of 
the phase (Appendix^, 



A: 



2 w^oo N 



1 ^ 



7j \'^3 



(20) 



J=0 



where F^^ (9) = 1 - 2h^^ (e) sin(26') + Ah^^ (e) cos^ 6 and 
7j stands for the species of the jth atom. Due to the 
periodicity of F{6) we can use the distribution functions 
of the phase to calculate the average H2U|) : 



A = 



2 ^ 

7,/3 



C7P7/3 



d0 



logFpi9)de 



- ^ log F^ (tt) - - ^ c^Pjfi / W, 

7 7,/9 



(21) 



From this last expression one can obtain numerically the 
Lyapunov exponent as a function of the energy. 
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III. RESULTS 

Once we have built the mathematical framework ly- 
ing under the model, one can try to answer the physical 
questions that obviously arise from this new configura- 
tion of the disordered quantum wire: how strong is the 
effect of this short-range correlations?, does it produce 
a measurable change in the density of states or in the 
localization of the electrons? As in Ref. I, we solve the 
functional equation numerically to obtain the density of 
states and to calcultate the Lyapunov exponent. It seems 
really hard to find analytical solutions for equations l|ll|) 
if they exist at all. Let us analize in the first place 
the density of states (DOS) for our quantum wires. In 
Figure 13 the evolution of the DOS for a binary wire as 
a function of the correlation pi2 is shown for a certain 
value of the concentrations. This distribution of states 
is drastically changed from the initial situation in which 
the probability to find the cluster -12- is low to the final 
stage when we impose the atoms of species 1 to appear 
always isolated. Note however that the concentrations 
are the same in both cases!. Another example can be 
seen in Figure 13 for a wire composed of 4 species. As 
the number of species increases, the number of correla- 





tions grows as 



m(m— 1) 



FIG. 2: Density of states for a disordered wire of two species 
with different correlations. ( — ) = —3 and ( — ) = 4 with 
ci = 0.4. 



and the distribution of states can 
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FIG. 3: DOS and Lyapunov exponent (grey line) for a 4 species quantum wire: 
tions are the same for the two configurations 



1, —1, 3, —5. Notice that the concentra- 




(a) ci = 0.4 



(b) P12 = 0.1 



FIG. 4: Variation of the Lyapunov exponent for a binary wire composed of [■^) = 1 and (^) 
correlation and as a function of the concentration. 



= —3, as a function of the 



adopt a lot of different shapes: the exploration of the 
whole correlations and concentrations space for m > 2 
can take a long time. From these graphics we conclude 
that the correlations can unbalance the spectra of the 
disorder system quite far from the completely random 
configuration. In fact the correlation can be tuned to 
open (close) an energy gap or to increase (decrease) the 
number of available states in a certain energy interval 
without changing the concentrations of the atomic com- 
ponents. The variation of the DOS with the correlations 
is of course a function of the concentrations being the 
chain with homogenized participation the one whose dis- 
tribution has the strongest dependence on P. Let us 
also remark that the fractal behaviour of the density of 
states in certain energy ranges, reported in reference I, 
still manifests itself for the different correlations regimes. 

Let us have a look at the localization of the electronic 
states. In Figures 0] and [3 the behaviour of the Lya- 
punov exponent for a binary chain changing the correla- 



tions and the concentrations is shown. In Figure 13 the 
different shapes of the degree of localization for different 
correlations can also be noticed for a 4-specics wire. The 
effective influence of the correlations on the localization 
of the electrons inside the quantum wires is established. 
In the study of correlated disorder two main kinds of 
correlations have been taken into account namely short- 
range and long-range. The former, mainly represented 
by the random dimer model^, has been widely accepted 
to be the reason for the emergence of resonant energies 
and a set of states close to the resonant one for which 
the localization length becomes larger that the system 
size4, improving greatly the transport properties^. The 
long-range correlations seem to be able to include in the 
spectrum mobility edges^, that can be entirely controlled 
by the correlators^, yielding an energy interval of a com- 
plete transparency. In our type of wires we have not 
observed mobility edges nor resonant energies. However 
an important effect on the localization of the states is 
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(a) ci = 0.5 



FIG. 5: Variation of the Lyapunov exponent for a binary w 
correlation and as a function of the concentration. 

shown that seems to act globaly in the whole energy 
range in contrast to the more restricted effect of others 
short-range correlated disorder models. In fact for cer- 
tain energies the Lyapunov exponent can be decreased 
an 80-90% of its maximun value changing the correla- 
tions at the expense of an increasing behaviour for other 
energies. Whether these correlations might cause A to go 
below the inverse of the length of a finite sample for some 
energy intervals is not certainly clear but it could be a 
point for future discussions. One can now argue if this 
influence of the correlations on the localization as well as 
on the DOS will be a measurable effect on short finite 
wires that represent a simple practical implementation, 
or on the contrary is just theoretically achievable only in 
infinite unreal systems. In order to answer that question 
and ascertain ourselves of the importance of the results 
we learn to generate finite sequences showing this kind of 
correlations. In FigureElwe show the Lyapunov exponent 
of an infinite wire for a certain concentration and corre- 
lation calculated using the functional equation, and the 
one obtained numerically from equation (|19|l for a corre- 
lated disordered chain of 1000 atoms. The convergence 
is excellent with only one realization of the correlated 
disorder using one thousand of atoms, and we have also 
checked it for others couplings of the species and other 
correlations obtaining a high degree of agreement. Thus 
it means that the described results are real effects with 
some practical importance on experimental realizations 
and on the other hand that the method of the functional 
equation, although designed for an infinite wire, can re- 
produce quite well the tendency in finite samples. 



IV. CONCLUDING REMARKS 

To summarize, we have extended the model proposed 
in reference I and completed the formalism introduced 
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(b) pi2 = 0.5 



:, composed of (^) = —2 and (^) = 5, as a function of the 
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FIG. 6: Lyapunov exponent for a binary wire with (^) = —2, 

(^) — 3, ci = 0.5 and pi2 — 0.9. The red line corresponds 
to the infinite wire. The blue line shown A numerically cal- 
culated for a sequence of 1000 atoms. 



there including the possible correlations which must be 
present in any realistic model of one dimensional conduc- 
tance. The method described here is the most natural 
manner to account for correlations that can manifest in 
nature or even inside a manufactured quantum wire in a 
nonintentional way. 

It has been established how due to the correlated dis- 
order the density of states is essentially modified and the 
localization of the electronic states changes in a non neg- 
ligible proportion which might alter some macroscopic 
properties of these structures. 

In addition to the obtained results we are preparing a 
future report in which we will construct the functional 
equation formalism in a more coherent way relating it to 
others methods currently used. There it will be shown 
how this mathematical framework is fiexible enough to be 
suitable for applying it to different disorder models such 
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as the Tight-Binding model as well as other correlation 
schemes such as the random dimer. 



If now we introduce the polar coordinates into the defi- 
nition of the Lyapunov exponent we obtain: 
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APPENDIX A: LYAPUNOV EXPONENT 

It can be shown that the electronic wave function sat- 
isfies, 



(Al) 



where Vfj is the amplitude of the wave function at the 
jth atomic site of the chain. This equation can be vi- 
sualized as a two dimensional mapping in polar coordi- 
nates with the identification ^'j+i = 7'j+icos0j+i and 



rj+i cos 9j+i 



2hy^{e) — 1\ frj cos 9 j 



1 







r ," smf 



(A2) 



which yields the known transmission relationship for the 
phase and the following for the moduli r, 



1 ^ 

A — lim — log 



J=0 



N 



— lim — log 



lim — y^ log 



N 



N^oo N 



3=0 



cos 9j+i 



3=0 



because according to the meaning of A the cosinus term 
must be zero. And using (|A3|) we can calculate the Lya- 
punov coefficient as a function of the phase, as: 



1 1 ^ 

A = - lim — Vlog^;.( 

2 N-^oo N ^ ^' 



(A5) 



= F^^ {9j) = 1 - 2h^^ (e) sm{29j) + Ah^ (e) cos^ 9j. 



(A3) which is our expression H2U|I . 
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